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摘 要 : 利用 精密 星 历 解 算 卫星 坐标 是 实现 GPS 高 精度 差分 定位 的 首要 基础 ， 
SIRF GPS 精密 星 历 的 采样 时 间 间 隔 ， 为 了 得 到 任意 时 刻 的 卫星 坐标 ， 需 要 对 
全 球 定位 系统 精密 星 历数 据 做 内 挂 和 外 推 处 理 。 分 析 卫 星 位 置 随时 间 变 化 的 规 
律 后 ， 采 用 广义 延 拓 法 对 卫星 的 一 段 轨道 构造 以 时 间 为 坐标 的 内 桂 和 外 推 模型 ， 
在 不 引入 较 大 位 置 误差 的 同时 简化 了 卫星 坐标 的 计算 过 程 ， 并 与 拉 格 朗 日 插值 
法 进行 对 比 实验 。 研 究 结 果 表 明 ， 利 用 广义 延 拓 法 内 醚 时 误差 小 于 5cm， 且 在 
外 推 30min 内 仍然 能 维持 精度 ， 明 显 优 于 拉 格 朗 上 日 插值 法 。 
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国际 GNSS 服务 中 心 (International GNSS Service, IGS) 提供 的 精密 星 
历 通常 用 于 全 球 定位 系统 精密 单 点 定位 等 数据 处 理 中 ， 但 国际 GNSS 服务 中 心 发 
布 的 精密 星 历 采样 间隔 是 15 min， 而 在 全 球 定位 系统 精密 定位 中 接收 机 的 采样 
率 一 般 远 小 于 这 个 值 ， 因 此 ， 需要 对 精密 星 历 进行 高 精度 的 插值 。 另 一 方面 ， 
由 于 国际 GNSS 服务 中 心 精密 星 历 仅 提 供 当 天 00: 00: 00°23: 45: 00 时 间 段 的 
星 历数 据 ， 若 想 仅 用 当天 的 数据 得 到 23: 45: 00°24: 00: 00 时 间 段 的 星 历 数 
据 ， 还 需要 进行 外 推 以 获得 任意 时 刻 的 卫星 坐标 。 

要 使 插值 定位 结果 和 平滑、 稳定 ， 要 求 插值 多 项 式 及 其 导数 连续 和 平滑 ， 较 
为 常见 的 插值 方法 有 Chebyshev 多 项 式 插值 法 、Newton 多 项 式 插值 法 和 
Lagrange 插值 法 等 “了”。 在 目前 的 数据 处 理 中 ，Lagrange 插值 应 用 较 广 ， 文 [4] 
指出 Lagrange 插值 法 在 本 质 上 与 Neville 插值 法 相同 ; 文 [5] 将 Lagrange 插值 
法 改进 ， 利 用 滑动 式 Lagrange 插值 法 获得 全 球 定 位 系统 精密 星 历 ; 文 [6] 对 比 
了 Lagrange 和 Newton 插值 法 的 内 插 效 果 ， 对 插值 结果 进行 了 总 结 分 析 。 上 述 
研究 中 ， 多 种 插值 方法 在 一 定时 段 内 需要 进行 高 低 阶 组 合 才 能 达到 最 优 效果 ， 
且 外 推 效 果 普 遍 不 很 理想 。 由 于 广义 延 拓 至 近 法 构造 的 插值 函数 在 所 有 同 阶 插 
值 函数 中 具有 最 小 平方 逼近 误差 ， 为 此 ， 本 文 将 基于 广义 延 拓 原理 对 全 球 定 位 
系统 精密 星 历 进行 插值 外 推 ， 利 用 求 得 的 逼近 函数 解 算 任意 时 刻 的 全 球 定位 系 
统 卫星 的 精密 坐标 。 


1 广义 延 拓 内 插 及 外 推 模 型 


广义 延 拓 逼近 法 在 分 片 边 表 点 上 满足 插值 条 件 ， 使 分 片 间 变化 协调 ， 并 充 
分 利用 分 片 插 值 区 域 的 周围 结 点 信息 ， 实 现 分 片区 域内 部 的 最 佳 拟 合 ， 有 机 结 
合 插值 法 和 拟 合法 的 优点 ， 在 数据 处 理 时 实现 了 高 精度 的 逼近 ”” 。 

首先 在 整 域内 进行 剖 分 处 理 ， 以 便 在 单元 域内 寻找 逼近 函数 及 基 函 数 ， 将 
定义 域 剖 分 成 个 互 不 重叠 的 子 域 : 

(1) 
单元 域 有 个 结 点 ， 对 应 结 点 坐标 为 ， 结 点 坐标 满足 定义 域 。 现 将 单元 域 延 


拓 至 ， 成 为 延 拓 域 ， 结 点 坐标 ,满足 延 拓 域 为 ， 如 图 1。 
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图 1 延 拓 域 及 相应 函数 值 
Fig. 1 Extension domains and their corresponding function values 


设 延 折 域内 有 个 节点 《〈《) ， 利 用 含有 较 多 结 点 的 延 拓 域 构造 单元 域内 逼近 
函数 : 
(2) 
其 中 ， 为 上 的 一 组 基 函 数 ; 为 待定 系数 ; 为 逼近 函数 项 数 ， 且 。 在 求解 函数 的 
待定 系数 时 ， 将 的 适用 范围 由 单元 域 扩 大 至 延 拓 域 。 则 在 延 拓 域 建立 如 (3) 式 
的 广义 延 拓 逼近 内 插 模型 ， 为 先 验 信息 采样 点 ， 为 逼近 函数 与 先 验 值 的 误差 : 
: (3) 
广义 延 拓 逼近 方法 将 单元 域 府 套 在 延 折 域 中 ， 吸 收 单 元 域外 延 拓 域 中 的 结 
点 信息 构造 单元 域内 的 拟 合 逼 近 函 数 ， 同 时 约束 单元 域 边 珊 结 点 ， 实 现 相 邻 单 
元 逼近 函数 的 协调 和 连续 ， 从 而 达到 单元 域内 还 近 函 数 的 最 佳 拟 合 效 果 。 利 用 
广义 延 拓 构造 的 逼近 函数 即 可 对 目标 变量 进行 插值 处 理 。 
广义 延 拓 插值 模型 不 但 可 以 用 作 内 插 模 型 ， 也 可 以 作为 外 推 模型 用 于 外 延 
应 用 的 场合 。 对 一 组 不 断 增长 的 数据 序列 ，， 已 知 以 前 的 数据 值 ， 根 据 先 验 数 
据 的 变化 规律 和 趋势 ， 欲 求 时 刻 的 值 ， 数 据 外 推 示 意图 如 图 2。 
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图 2 数据 外 推 示意 图 


Fig. 2 Data extrapolation intent 
按照 广义 延 拓 插 值 外 推 的 设计 理念 ， 令 为 最 新 时 刻 ， 采 用 外 推算 法 可 得 下 
一 时 刻 的 值 。 建 立 广义 延 拓 外 推 模 型 : 


(4) 
在 上 述 模型 中 ， 由 于 插值 点 为 最 新 采样 点 的 值 ， 为 了 克服 这 一 点 数值 突变 
带 来 的 误差 ， 考 虑 建立 改进 的 广义 延 拓 模型 : 
. (5) 
即 选择 最 新 的 个 采样 点 的 平均 值 作为 约束 ， 改 进 模型 的 解法 与 原 模型 的 解 
法 基本 一 致 ， 但 数据 的 平稳 性 更 好 ， 插 值 拟 合 曲线 更 平滑 ， 外 推 精度 更 高 。 


2 全球 定位 系统 精密 星 历 的 内 插 和 外 推 方法 


2.1 全 球 定位 系统 精密 星 历 的 解 算 方法 

考虑 卫星 的 实际 运行 轨道 会 偏离 无 摄 运行 轨道 ， 因 此 全 球 定位 系统 的 星 历 
参数 通常 采用 一 套 扩展 后 的 开 普 勒 轨道 参数 ， 由 16 个 参数 组 成 ， 包 括 6 个 开 普 
勒 根 数 ，9 个 摄 动 参数 和 星 历 参考 时 间 。 

由 文献 [10] 提 供 的 解 算 算法 ， 可 将 卫星 定位 计算 过 程 简要 描述 如 下 : 


， (6) 
其 中 ，，， 为 每 个 历 元 待 求 的 卫星 坐标 。 经 过 坐标 变换 和 摄 动 校正 ， 最 终 得 到 
t 时 刻 该 卫星 在 WGS-84 地 心地 固 坐 标 系 中 的 坐标 计算 矩阵 : 
(7) 


其 中 ， 为 信号 发 射 时 刻 卫星 在 轨道 平面 的 位 置 坐标 ; 为 信号 发 射 时 刻 的 升 交点 
赤 经 ; 为 轨道 倾角 ， 均 可 由 星 历 参数 计算 求 得 。 

由 以 上 求解 过 程 可 知 ， 利 用 星 历 参 数 计算 全 球 定位 系统 卫星 在 某 一 时 刻 的 
空间 位 置 必 须 提供 准确 的 星 历数 据 ， 每 定位 一 次 需要 的 计算 量 也 挺 大 。 
2. 2 全 球 定位 系统 精密 星 历 的 插值 可 行 性 分 析 

图 3 中 的 三 条 曲线 分 别 是 某 颗 卫星 的 位 置 在 WGS-84 地 心地 固 坐标 系 中 的 
义 Y 和 2 分量 随时 间 的 变化 情况 ， 历 元 间隔 为 15min， 可 见 卫星 的 轨道 位 置 呈 周 
期 性 变化 。 如 图 4， 卫星 位 置 的 各 分 量 在 一 段 短 时 间 内 随时 间 变 化 平滑 ， 几 乎 
呈 线 性 变化 。 
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图 3 长 时 间 内 卫星 的 空间 位 置 图 4 短 时 间 内 卫星 的 空间 位 置 
Fig. 3 Space position of satellite over a long period of time Fig. 4 Space position of 


satellite in short time 


因此 ， 了 卫星 的 一 段 轨道 可 以 用 一 个 以 时 间 域 上 的 插值 多 项 式 表 达 "" 1， 如 果 
此 逼近 函数 构造 得 当 ， 那 么 这 种 插值 方法 不 会 引入 很 大 的 卫星 位 置 误差 ， 可 避 
免 2.1 节 所 述 的 复杂 计算 方法 ， 使 计算 量 大 幅度 减少 。 

2. 3 基于 广义 延 拓 的 全 球 定 位 系统 精密 星 历 内 插 和 外 推 方法 

利用 广义 延 拓 插 值 原理 ， 对 卫星 位 置 与 速度 在 WGS-84 地 心地 固 坐标 系 中 的 
X, Y, 2 分 量 分 别 分 段 建立 广义 延 拓 模型 ， 在 到 这 段 时 间 内 ， 以 卫星 位 置 在 X 方 
向 上 的 分 量 为 例 ， 若 在 区 间 内 取 两 端点 作为 先 验 点 进行 约束 ， 则 广义 延 拓 插值 
模型 为 

; (8) 
其 中 ， 表 示 在 插值 区 间 左 端点 时 刻 卫星 位 置 在 WGS-84 地 心地 固 坐标 系 中 XY 方 向 
上 的 分 量 ; 表示 在 插值 区 间 右 端点 时 刻 卫 星 位 置 在 WGS-84 地 心地 固 坐 标 系 中 XX 
方向 上 的 分 量 ; 为 待 求 系数 ; 是 精密 星 历 更 新 的 时 间 点 。 


对 上 述 模型 展开 求解 ， 即 由 
(9) 
得 到 法 方程 ，〈9) 式 中 ， 为 拉 格 朗 日 乘 子 ， 其 中 


A, ， 解 上 述 法 方程 ， 求 得 系数 后 ， 可 得 卫星 位 置 在 X 方 向 上 的 分 量 随时 
间 变化 的 广 义 延 拓 内 插 函 数 ， 将 所 需 时 刻 时 间 代 入 ， 即 可 快速 求解 。 
接着 ， 由 (5S) 式 构造 广义 延 拓 外 推 模型 : 


41) 
其 中 ， 变 量 含义 见 8) 式 ， 其 中 表示 选用 了 最 新 个 卫星 位 置 的 《分 量 ， 取 平均 
人 对 此 模型 的 解法 如 下 : 


(412) 
可 得 系数 求解 矩阵 ; 
(13) 
Hh: 
(14) 


解 得 待 求 系数 后 ， 便 可 得 到 外 推 公式 : 


(15) 
同 理 可 得 卫星 位 置 在 其 余 各 分 量 上 的 外 推 公式 。 


3 算 例 分 析 


现 验 证 广义 延 拓 插值 外 推 模型 的 推算 效果 ， 采 用 IGS 提供 的 2017 年 6 月 6 
日 ( 即 GPS 1952 周 ) 采样 时 间 间 隔 为 15min 的 精密 星 历 数据 ， 选 用 的 卫星 为 
PRN15 号 。 选 取 前 4 个 采样 数据 作为 移 验 采样。 点 ， 以 Lagrange 插值 法 作为 参照 ， 
并 与 精密 星 历 数据 拟 合 的 标准 曲线 对 比 。 得 到 两 种 精 密 星 历 插值 法 的 逼近 效果 
如 图 5。 
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图 5 两 种 精密 星 历 插值 法 的 插值 外 推 效 果 对 比 


Fig. 5 Comparison of interpolation and extrapolation effects of two precise ephemeris interpolation 
methods 


由 图 5 可 见 ， 当 以 1 小 时 的 精密 星 历 作为 先 验 点 进行 拟 合 逼 近 函 数 时 ， 在 
1.5 小 时 的 时 段 内 ，Lagrange 插值 法 和 广义 延 拓 法 构造 的 逼近 函数 与 精密 星 历 
的 拟 合 曲线 基本 吻合 ， 而 在 1.5 小 时 以 后 ， 广 义 延 拓 法 的 外 推 效果 明显 好 于 
Lagrange 插值 法 。 为 了 更 直观 和 导 观 察 广义 延 拓 法 构造 的 模型 精度 ， 将 X、Y，2Z 
各 分 量 的 误差 整合 为 位 置 误差 6。，， 得 到 位 置 误差 随时 间 变 化 如 图 6。 
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图 6 卫星 位 置 随时 间 变 化 图 


Fig. 6 Chart of satellite position versus time 


由 图 6 可知， 前 4 个 比 对 点 拟 合 的 曲线 为 广义 延 拓 内 插 效 果 ， 其 插值 的 卫星 
位 置 误差 小 于 5cm， 自 第 5 个 比 对 点 开始 为 广义 延 拓 的 外 推 效 果 ， 可 见 广 义 延 拓 
外 推 1 小 时 左右 的 卫星 位 置 误差 小 于 10 cm， 外 推 2 小 时 的 卫星 位 置 误差 小 于 
20cm， 而 在 外 推 2 小 时 以 后 ， 卫 星 位 置 误差 开始 迅速 变 

由 于 16S 提供 的 精密 星 历 误差 精度 为 5cm， 因 此 ， 利 用 广义 延 拓 法 对 全 球 定 
位 系统 精密 星 历 的 内 插 能 够 满足 精度 要 求 ， 而 在 利用 广义 延 拓 外 推 时 ， 在 30 
min 内 也 能 维持 精度 ， 在 较 长 时 间 内 也 不 会 引入 很 大 的 位 置 误差 。 


4 结论 


采用 广义 延 拓 法 对 全 球 定位 系统 精密 星 历 进行 卫星 位 置 内 插 时 效果 较 好 ， 
可 以 将 卫星 的 一 段 轨道 用 一 个 以 时 间 为 坐标 的 模型 表达 ， 使 得 计算 卫星 位 置 时 
在 不 损失 精度 的 情况 下 ， 不 再 进行 复杂 元 长 的 步 又 ， 极 大 程度 减少 了 计算 量 。 
此 外 ， 通 过 实例 分 析 ， 广 义 延 拓 外 推 法 能 够 保证 外 推 1 小 时 的 精度 能 满足 要 求 ， 
广义 延 拓 内 插 法 可 获得 全 球 定 位 系统 任 一 时 刻 的 精密 星 历数 据 。 
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Application of Generalized Continuation Method in 


Interpolation and Extrapolation of GPS Precise Ephemeris 
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Abstract: The use of precise ephemeris to solve satellite coordinates is the primary basis for 
high precision differential positioning of GPS, in order to obtain satellite coordinates at any time 
since the GPS precise ephemeris has limit of the sampling interval, it is necessary to interpolate 
and extrapolate GPS precise ephemeris data. After analyzing the law of satellite position changing 
with time, the interpolation and extrapolation model which is relat to time is built by generalized 
extension method, and it is used to describe a section of satellite orbit, the calculation process of 
satellite coordinate is simplified without introducing larger position error through this method, 
then compared with the Lagrange interpolation method in experiment. The result shows that the 
error of interpolate when using generalized extension method is less than 5cm, and maintain the 
accuracy of extrapolation in 30 minutes, significantly better than the Lagrange interpolation 
method. 
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